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Context. Limb darkening is an important tool for understanding stellar atmospheres, but most observations measuring limb darkening 
assume various parameterizations that yield no significant information about the structure of stellar atmospheres. 
Aims. We use a specific limb-darkening relation to study how the best-fit coefficients relate to fundamental stellar parameters from 
spherically symmetric model stellar atmospheres. 

Methods. Using a grid of spherically symmetric Atlas model atmospheres, we compute limb-darkening coefficients, and develop a 
novel method to predict fundamental stellar parameters. 

Results. We find our proposed method predicts the mass of stellar atmosphere models given only the radius and limb-darkening 
coefficients, suggesting that microlensing, interferometric, transit and eclipse observations can constrain stellar masses. 
Conclusions. This novel method demonstrates that limb-darkening parameterizations contain important information about the struc- 
ture of stellar atmospheres, with the potential to be a valuable tool for measuring stellar masses. 
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1. Introduction 

Limb darkening, the change of surface intensity from the center 
to the edge of a stellar disk, is a powerful measure of the physi- 
cal structure of stellar atmospheres. However, it is difficult to ob- 
serve the intensity profile across a stellar disk except for the Sun. 
Most observations of stellar limb darkeni ng are indirec t, coming 
from light cur ves of eclipsing bin aries ((Claret 2008) or plan- 
Knutson et al. I f2007h . interferometric visibilities 
2009) ng light curves dAn et al.l 



l CUIV 

isits f] 
et all 



etary nansits 
(Haubois 

2002), and these indirect me thods have limited precision (e.g. 
Poppeil l 19841 IZub et al.ll20TTl) . Because of this, limb darkening 
is commonly treated as a parameterized function of fi = cos 8, 
where 9 is the angle between the direction to the distant observer 
and the normal direction at each location o n the stellar surface . 
An example of a linear parameterization is dSchwarzschildll 906) 



hill = 1) 



1 -a A (l -/j). 



(1) 



As numerical model atmospheres became robust, the com- 
puted intensity profiles were represented by more elaborate 
limb-darkening parameterizations that include higher order 
terms of u or are expressed in terms of powers of r — sin# 
dHevrovskvll2007l) . as well as being normaliz ed with respect to 
the st ellar flux instead of the central intensity dWade & Rucinskil 
119851) . These more detailed limb-darkening laws were used to 
interpret the improving observations noted above. 

However, even with these advances, current limb-darkening 
observations are still unable to constrain the model stellar 
atmospheres. Recent interferometric observations of nearby 
red giants are not yet precise enough to differentiate be- 
tween the center-to-limb intensity profiles from Phoenix and 
Atlas model atmospheres or even between plane-parallel or 



spherically symmetric m odels (Wittkowski et al. 2004, 2006a b; 
iNeilson & Lesteti 1200 8!). However, the combination of the ob- 
servations and models do provide constraints for fundamental 
stellar parameters such as the effective temperature and gravity. 

In this work, we continue our previous analysis 
(INeilson & Lester! 1201 ll hereafter Paper 1) to explore the 
connection between limb darkening and stellar parameters. We 
adopt the flux-conserving law 



h(M) 



1-A A 1 



where "H\ is the Eddington flux defined as 



'J 



(2) 



(3) 



because this law was used bv lFields et al. d2003h to ana lyze the 
microlensing event EROS-BLG-2000-5 (|AiLetd]|200l- We fit 
this limb-darkening law to the intensities of spherically symmet- 
ric model stellar atmosp heres computed with the SAtlas code 
(Lester & Neilson 2008) to explore the relation between funda- 
mental parameters and limb-darkening coefficients hinted at in 
Paper 1. 



2. Limb darkening and fundamental parameters 

Many previous studies have used the limb darkening predicted 
by model stellar atmospheres as a tool to achieve a better deter- 
mination of a stellar diameter dMerand et all [2 10) or a better 
characterization of a transiting planet dLee etalJ2 012). Our goal 
is different; we want to determine how limb darkening, repre- 
sented by Eq. (0, depends on the fundamental parameters of 
stellar atmospheres. Because there are several steps leading to 
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Fig. 1. Best-fit V-band limb-darkening relations using Eq. (|2]i for 
grids of 5 Mq spherically symmetric and plane-parallel Atlas 
model atmospheres. A fixed /u-point is clearly seen in both cases. 



our conclusion, we describe our approach step-by-step to make 
it as clear as possible. 

2. 1 . Intensity fixed point 

Figure [TJ plots the curves produced by the limb-darkening law 
given by Eq. (|2]i for the V-band intensities computed for both 
plane-parallel and spherical Atlas models. The cube of spherical 
models is from Paper 1, with luminosities and radii correspond- 
ing to the range of 3000 K < T eS < 8000 K in steps of 100 K, 
-1 < logg < 3 in steps of 0.25 and 2.5 < M/Mq < 10 in steps 
of 2.5 Mq, which we have supplemented with additional mod- 
els having masses of M — 0.5 and 1.0 Mq. The grid of plane- 
parallel models spans the equivalent range of effective temper- 
ature and gravity. Figure [TJ shows that both sets of models have 
a fixed point, fit, where the curves intersect, although there is a 
larger spread for the spherical models because of the combina- 
tion of the three fundamental stellar parameters, L*, M* and 
The immediate vicinity of the spherical fixed point is displayed 
in more detail in Fig. 6 of Paper 1. 

2.2. Relationship between A and B 

In Paper 1 we showed that the fixed point, n\, is caused by the 
limb-darkening law's Ax -coefficient being linearly related to the 
Bx -coefficient, Ax °c cc\Bx- The proportionality term, ax, de- 
pends only on two quantities, both of which are functions of the 
stellar intensity: the mean intensity, Jx, and the pseudo-moment 
fx, defined as 



Combining these two radiation terms defines a new quantity, 



which is similar to the Eddington factor, f — K/J dMihalasI 
[l978h . As we showed in Eq. (19) of Paper 1, ax depends only 
on 7]x, 

(5^/4- l)/6+ 11/96 
" A (5^/4 -l)/4+ 1/6 ' 

2.3. Relationship between a and n\ 

Our analysis in Paper 1 also showed that the fixed point, is a 
function of ax- Rearranging our result from Paper 1, we obtain 

3a A (l + ax) + f| + au(l +a A ) + §§ 
V\ = — 2 ■ (7) 

Working back through this chain of logic, we find that because 
ji\ depends on ax, which depends on rjx, we are led to the conclu- 
sion that Ix(Mi)/ ^Ha depends solely on the two angular moments 
of the stellar intensity, Jx and fx- 

The analytic derivations of Eq. (|6]l and Eq. (0 in Paper 1 
assumed, for convenience, that r\x is the same for all models. 
This assumption is true for plane-parallel models because they 
share the same flat geometry. However, r\x varies for spherical 
atmospheres because they have different amounts of curvature 
depending on the values of their particular fundamental param- 
eters. As a result, the values of ax and \i\ also vary for spherical 
atmospheres. However, this variation does not alter the logical 
connections between these quantities, as we show numerically 
in the following section. 

2.4. Dependence on surface gravity 

Having used our analytic results to trace the logical connec- 
tion between the key variables, we now continue our exploration 
using large grids of plane-parallel and spherical model atmo- 
spheres to compute the key quantities numerically, which en- 
ables us to drop the assumption that rjx is constant. Figure [2] 
shows the results for the V-band. 

Beginning with the plane-parallel models, the black crosses 
in the upper left plot of Fig. [2] shows that the 77 parameter is 
nearly constant over a very wide range of surface gravity for 
the models being considered. As required by Eq. (0, the near 
constancy of 77 causes the a coefficient also to be almost constant, 
as shown by the black crosses in lower left panel of Fig. [2] Next, 
as a result of Eq. (|7), the fixed point, flu must be nearly constant, 
as shown by the black crosses in the upper right panel of Fig. [2] 
Finally, when the nearly constant fi\ is used in Eq. ©, the result 
is a nearly constant value of the normalized intensity at the fixed 
point. This conclusion is confirmed by the black crosses in the 
lower right plot of Fig. [2] 

The spherically symmetric models are fundamentally differ- 
ent. The various colored symbols in the upper left plot in Fig. [2] 
confirms the result of Sect. 12.31 that the parameter 77 does vary 
with surface gravity. It then follows from Eq. © that a must 
also vary, as shown by the colored symbols in the lower left plot. 
As before, Eq. (0 requires that the variation of a leads to the 
variation of [i\ shown by the colored symbols in the upper right 
plot of Fig. [2] Finally, this leads to the variation of the normal- 
ized intensity at the fixed point shown by the colored symbols 
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Fig. 2. Surface gravity dependence of r\x, ax, fixed point fi\ and the normalized intensity at the fixed point, I\(n\)l1'H\, for the V 
spectral band. Plane-parallel model atmospheres are denoted by black crosses, and spherical model atmospheres are represented 
by other symbols. Filled red circles are T e s = 3000 K, open green squares 4000 K, open blue circles 5000 K, downward pointing 
magenta triangles 6000 K and upward pointing pale blue triangles are 7000 K. 



in the lower right plot. For spherically symmetric models, the 
conclusion is that the location and intensity of the fixed point 
depend on the fundamental stellar parameters. 

The beginning step of this logical chain, the dependence of 
T] on the atmospheric parameters, is physically reasonable. As 
shown in the upper left panel of Fig. [2] as logg decreases, 77 in- 
creases towards unity, meaning that P —> J. This is the result of 
the intensity profile becoming extremely centrally concentrated, 
or, equivalently, the atmospheric extension approaches the stellar 
radius. Because rj depends upon the extension of the atmosphere, 
so do a, the fixed point fii, and the intensity at the fixed point. 



3. Parameterization of atmospheric extension 

To explore further the results shown in Fig. [2] we replace \ogg 
by a more explicit representation of the extension of the stellar 
atmosphere, for which we adopt 



(8) 



This dependence is physically plausible because the extension 
is larger for stars of larger radius, but greater masses pull the 
atmospheres into a more compact configuration. 



One way of deriving Eq. (O begins with the definition of the 
pressure scale height, 



H = 



kJ 
fim H g' 



(9) 



where H is the distance over which the pressure changes by a 
factor of e, k is the Boltzmann constant, T is the gas tempera- 
ture, fj. is the mean molecular weight of the gas, mu is the mass 
of the hydrogen atom, and the surface gravity g = GM+/R+. 
Using H to represent the thickness of the atmosphere, the rela- 
tive extension of the atmosphere is 



AR _ H 
~R^~~R+ 



(10) 



both sides of which are dimensionless. Assuming that T and fi 
are nearly constant over the distance H, the relative extension of 
the atmosphere becomes 



AR 
^7 



constant 
R+MJRl 



(ID 



which is our expression. 

To test the correlation of R+/M+ with the actual atmospheric 
extension, we use the Rosseland optical depth to determine 
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Fig. 3. Extension parameter /?*/M+ compared to the relative at- 
mospheric extension, AR/R+, derived for more than 3000 spher- 
ical models. 

for each of our several thousand spherical models the quantity 
AR/R+, where AR is the physical distance from the location 
where tr oss = 1 out to the location of where tr oss = 0.001, which 
is slightly more than one pressure scale height, and is the ra- 
dius where tr oss = 2/3. Figure [3] plots the extension parameter 
R+/M+ for each model versus that model's AR/R+. While there 
is a spread because T and yU do vary over the distance H, it is ob- 
vious that R+/M+ tracks the actual extension of the models quite 
well. 

Using R+/M+ to represent the atmospheric extension trans- 
forms Fig. |2] into Fig. 21 where the extension parameter for 
the plane-parallel models is determined by assuming a mass of 
5 Mq and computing the radius for each model from its surface 
gravity. Figure [4] shows, as expected, that r\x is a tight function 
of log(/?*/M*), except for the models with T eB = 3000 K. It 
follows, because of the connections established in Sec. |2l that 
all other quantities also depend on R^/M*. For plane-parallel at- 
mospheres the quantities are nearly constant, as expected, and 
as they were found to be in Fig. [2] Therefore, the values of //i 
and Ia(m\)/2'Ha depend on the fundamental stellar parameters 
R+ and M*. The deviation of the models with T e g = 3000 K 
may be due to changes in the atmospheres at the lowest temper- 
ature. 

Excluding the models with r e ff = 3000 K, we fit for each 
waveband the dependence of n\ and I / \{fii)/2'J-( / \ on log(/?*/M*) 
shown in Fig. [4] using the functional forms 



Pi=C, log^ 



and 



+ A 



(12) 



(13) 



The values of C M , D^, C\, and D\ for these fits to the spherical 
models are given in Table Q] for the five wavebands B, V, R, I, 
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Fig. 5. V-band values of the A limb-darkening coefficient (top 
panel) and the B limb-darkening coefficient (bottom panel) for 
the linear-plus-square-root parametrization given by Eq. (0, 
plotted as a function of the extension parameter, R+/M+ in solar 
units. The symbols have the same meaning as in Fig. [2] 



and H. We include fits for models with effective temperatures 
r e jf = 5000 K and models within the effective temperature range 
of 4100 - 4300 K (r eff = 4200 + 100 K). Equation CES) and 
Eq. (fT3l l will be employed in later sections. 



4. Iterative method 

The right panels of Fig. |4]show that fii and I / i(jj[)/2'H / i depend 
on R+/M+ for the spherical model atmospheres. It follows from 
the limb-darkening law, Eq. (0, that the coefficients and B\ 
must also be functions of the atmospheric extension, which is 
confirmed in Fig. [5] Because the limb-darkening coefficients 
and B \ in Eq. (0 can be derived from observations, this estab- 
lishes an observational method for determining the atmosphere's 
extension. 

In Fig. |5] the open colored symbols show that the limb- 
darkening coefficients of stars with > 3000 K share 
a common dependence on the atmospheric extension for 
log(7?*/M*) ^ 1.5 in solar units. This suggests it is possible 
to derive a unique value of R+/M+ for stars with atmospheric 
extensions R*/M + £ 30 Rq/Mq. Stars with log(^*/M*) < 1.5 
show a clear dependence on T e s, 
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Fig. 4. (Left) Values of the pseudo-moment rjy and variable ay as a function of atmospheric extension for spherical atmospheres, 
R-t/M* in solar units. (Right) The top panel shows the value of the primary fixed point, fi\, for the linear-plus-square-root 
parametrization, also plotted as a function of the extension of the spherical atmosphere. The bottom panel shows the dependence of 
the normalized V-band intensity of the fixed point on the atmospheric extension. The symbols have the same meaning as in Fig. [2] 
When comparing to Fig. [2] note that increasing extension corresponds to decreasing logg. 



The tight dependence of the yUi and Ia(ji\)/2 < Ha on 
displayed in Fig. [4] follows directly from Sect. I2.31 where we 
showed that /Ji is a function of ax, which depends only on 
r\x = Pa/ J a- The more scattered dependence of the coefficients 
A a and Ba on the extension shown in Fig. [3] is due to these co- 
efficients having a less direct functional dependence on the ex- 
tension through the separate terms Ja/2'Ha and 'PaI'^Ha, as we 
showed in Paper 1 . Because of this, as we will show in the fol- 
lowing sections, the values of Aa and Ba do not measure atmo- 
spheric extension as cleanly as using \i\ and Ia([1\)I2 < Ha as con- 
straints. 

Although Aa and Ba are less direct measures of the atmo- 
spheric extension, they are determined by the modeling of the 
lensing observations. Therefore, we propose an iterative method 
to use Aa and Ba to determine the atmospheric extension. The 
method begins by assuming an initial value for log(/?*/M*), 
which is used to compute \i\ with Eq. d!2t . This value of [i\ 
and the coefficients of the fit given in Table Q] enable us to com- 
pute the normalized intensity at the fixed point using Eq. (|2]). 
Finally, we compute a new value of log(/?*/M*) using Eq. ( fT3l ). 
This new value of log(7?*/M*) replaces our initial assumption, 
starting a second iteration, which is continued until the value 
of log(/?*/M*) converges. If the stellar radius is determined by 



other methods, such as from interferometric observations and 
distance estimates, the converged value of log(R*/M*) can be 
solved for the star's mass. 



5. Test using a model atmosphere 

The previous section describes a potential method for deter- 
mining the fundamental parameters of a star using stellar limb- 
darkening laws and the underlying physics behind these laws. 
To test this method, we arbitrarily selected a model with r e ff = 
5000 K, logg = 2, M* = 5 Mq, corresponding to ^*/M* = 
7.42 Rq/Mq or log(fl*/M*) = 0.87. Note that this value of 
the atmospheric extension falls in the region of Fig. where 
the limb-darkening coefficients also depend on the T e ff, whereas 
the fixed point and intensity in Fig. [4] are less dependent on 
the r e ff. By selecting a trial model for which the Aa- and B\- 
coefficients depend on as well as the atmospheric extension, 
we can determine how the iterative method is affected by this ad- 
ditional dependence. Table|2]lists the test model's A and B limb- 
darkening coefficients for each of the five wavebands. Using the 
iterative method outlined in Sect. [4] we computed the /?+/M* 
for each waveband, with the results given in Table [2] By assum- 
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Table 1. Best-fit coefficients for Eq. (TTZl l and Eq. ( fT3l l as a function of log(R*/M*) for various wavebands and effective temperatures. 
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fundamental parameters given earlier; this leads to the spread of 
points shown in Fig. [4] 

The variation with wavelength also contributes to the uncer- 
tainty. Averaging over all wavelengths, the mean difference of 
the predicted R^/M* from the model's value is only 5%, with 
the largest deviation being about 13% for the V-band. The wave- 
length variation of /?*/M* is related to the definition of the stel- 
lar radius, which we define as = ^(tr oss = 2/3). However, 
in Table [2] the different wavelength bands produce values of 
R+ - ^(tb,v,r,i,h = 2/3), which will differ slightly from the 
radius where tr oss =2/3. 

Using the five wavebands in Table [2] we find CR*/M*) avg = 
7.01 ± 0.45 Rq/Mq, and logg avg = 2.05 ± 0.05, both of which 
are consistent with the actual values. Therefore, our proposed 
method for using limb-darkening coefficients recovers the actual 
stellar parameters for the theoretical test case with a statistical 
uncertainty of w 8% in R+/M+. 

We improved the test by fitting simultaneously the limb- 
darkening coefficients for all the wavelength bands using the 



Table 2. Limb-darkening parameters and predicted R*/M* for 
a model stellar atmosphere with T e g = 5000 K, logg = 2 and 
M = 5 M Q . 



Band 


A 


B 


RJM+ 


l°g# 


B 


0.613 


0.507 


6.65 


2.09 


V 


0.048 


1.206 


6.46 


2.12 


R 


-0.300 


1.622 


7.10 


2.04 


I 


-0.642 


2.051 


7.28 


2.02 


H 


-1.189 


2.719 


7.54 


1.99 



ing a stellar radius of = 37.1 Rq, we determined the surface 
gravity from /?*/M*, which is also given in Table [2] 

Even though the values of and for our test model have 
no observational uncertainty, the results in Table|2]differ slightly 
from the test model's values. One reason for this difference is 
that the coefficients given in Table [TJ have uncertainties because 
they are fits to a grid of thousands of models with the range of 
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fix-point relations for all models with T e ff > 3000 K. This re- 
sulted in the best-fit values /?*/M* = 6.83 + 0.38 Rq/Mq, and 
logg = where the error bars come from the formal 

uncertainty of the fits to the coefficients in Table Q] which are 
used in Eq. dT2b and Eq. ( fT3l . When the fit was repeated using 
only the relations for models with T e ff = 5000 K, the result was 
R+IM+ = 7.6r°;| Rq/Mq and logg = We conclude 

from this test using a model stellar atmosphere that our proposed 
method for measuring stellar properties yields consistent results, 
although, as expected, the uncertainty of /?*/M* increases for 
stars with smaller atmospheric extension. 

6. Conclusions 

The purpose of this study has been to investigate how limb dark- 
ening probes fundamental stellar properties. The limb darkening 
was parameterized by a flux-conserving linear-plus-square-root 
law that has two free parameters that are functions of two angular 
moments of the intensity. The ratio of these two moments pro- 
vides a measure of the amount of stellar atmospheric extension, 
represented by R+/M+. We tested our method for determining 
the extension from limb-darkening parameters using a spheri- 
cally symmetric model atmosphere and found agreement of the 
mean value of the five spectral bands analyzed to within 8%. 

We also find significant variation of R+/M+ due to the def- 
inition of the stellar radius. From the context of the model of 
the stellar atmosphere, we find it natural to define the stellar 
radius where tr oss = 2/3, but other options are possible, such 
as Tr oss = 1 or using a specific wavelength, such as 500 nm, 
to define the radius where T500 = 1. Observations, of course, 
are done in specific wavebands, such as V, and the stellar ra- 
dius is assumed to refer to some particular optical depth, such as 
Ty = 2/3. Switching to another waveband will lead to another 
value of the radius. As an example, in the near-infrared the pre- 
dicted radius is larger than at optical wavelengths. Therefore, the 
definition of the stellar radius is ambiguous in general for stud - 
ies of angular diameters (e.g. I Wittkowski et alj r2004. 2006b a), 
which makes the definition of the stellar radius challenging for 
st ars with signif icant atmospheric extension. A recent example 
is Ohnaka et al.l d201 lh . who measured the /f-band angular di- 
ameter of B etelgeuse to be about 2 mas (« 5%) smaller than 
measur ed bv lHauboisetaLl d2009t) in the //-band. However, si- 
multaneously fitting R+/M+ to multiband limb-darkening obser- 
vations provides a bolometric value for /?* /M* analogous to how 
spectrophotometry can be used to measure a star's effective tem- 
perature. 

The connection between atmospheric extension and stel- 
lar limb-darkening is a consequence of the assumption of 
spherical symmetry for the stellar atmospheres, not a unique 
feature of our SAtlas code. The results presented here 
should also be found using intensity profiles from spherically- 
symmetric Marcs models and the Ph oeni x model atmosphere s 
used by IClaret & Hauschildtl (120031) and iFields et all (120031) . 
Unfortunately, in both of these works the authors truncated the 
intensity profiles to remove the low intensity limb and then re- 
distributed the spacing of /i-points. Truncating the intensity pro- 
file eliminates critical information about the limb darkening, and 
makes the star's intensity distribution resemble a plane-parallel 
atmosphere, as we demonstrate in Fig. [6] 

Our method for determining fundamental stellar parame- 
ters is based on the general atmospheric properties of flux- 
conserving limb-dar kening law s, not the pa rticular law used 
here. For instance, iHevrovskvl (120031) and IZub et all d2011l) 
found a fixed point in flux-conserving linear limb-darkening 



laws, which indicates that the linear limb-darkening coefficient 
is a measure of the mean intensity and flux of the atmosphere. 
Another example is a quadratic limb-darkening law where the 
limb-darkening coefficients would provide a measure of the 
mean intensity and second moment of the intensity, K. Because 
these laws measure the correlation between various moments of 
the intensity in an atmosphere, they would also measure the at- 
mospheric extension of that star. 

We conclude that the method outlined in this work is a viable 
way to use a limb-darkening law, such as the linear-plus-square- 
root law employed here, to determine the fundamental physi- 
cal parameters of stars. The method requires knowledge of one 
of the degenerate parameters effective temperature or luminos- 
ity. However, the method is able to constrain stellar parameters 
using limb-darkening observations even in the case when those 
limb-darkening observations from microlensing or eclipsing bi- 
naries cannot test the model stellar atmosphere directly. It seems 
as though we still have things to learn from fairly simple repre- 
sentations of limb darkening, and that, as observations continue 
to improve, this will become an even more powerful tool in the 
study of stellar astrophysics. 
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Fig. 6. (a) Surface-brightness distribution as a function of fractional radius. The solid line is the full set of intensities from a spherical 
model atmosphere with T s $ = 5000 K, logg = 2.0 and M* = 5 Mq. The dashed line has removed intensity values for r/R* > 0.995 
and then renormalized the fractional radius to the interval 0—1. The dotted line has removed intensity values for r/R+ > 0.98 and 
then renormalized the fractional radius to the interval - 1 . (b) Surface-brightness distributions from panel (a) plotted as a function 
of yu, with the lines having the same meaning as in the (a) panel, (c) Fits to the surface-brightness distributions in panel (b) using 
Eq. (0. The lines have the same meaning as the other two panels. 
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